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Elastic Stability and Deflections of Eccentrically - 
Loaded Compression Members Made of Material With No 
Tensile Strength 


by 


Felix Y. Yokel 


A mathematical solution is derived, which permits the 
computation of critical load, deflection and stresses for 
eccentrically loaded slender prismatic compression members, 
made of materials that develop compressive strength but have 
no tensile strength. A graphical presentation of the solu- 
tion facilitates its application. In an example of applica- 
tion, the solution is used to compute the strength of masonry 
walls which were tested by the Structural Clay Products 


Institute. There is good agreement between computed and mea- 
sured strength. 


Key Words: Buckling, compression members, deflection, 
equilibrium, load eccentricity, masonry, 


stability, section cracking, stress distribution, 
unreinforced concrete. 
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ie introductian and Scope 

itive paper SOluttons are déerived,for the computation 
of elastic deflections and stability of compression members 
made of materials with compressive strength and no tensile 
strength and with a linear relationship between compressive 
stresses and strains. Computed values could be considered 
Woe enn Sto det hectr ons, and lower limits for critical 
tieoad ot (sSolvd prismatie walls made of masonry or unrein- 
HOnGCC COMecLece. Since these materials have properties simi- 
ih sO elosceassumed im sitive paper. in Appendix Il1 the solu- 
PLoS sane sUSecd tOupcedicte the results of 39 full-scale tests 
on buick walls which were conducted by the Structural Clay 
Products Tages anes Toya and the correlation between test 


PeoctiLes iid theory ets discussed,. 


7a Assumed Material Properties and Member Geometry 

The material is assumed to be Biasizie, Go ynave a lanear 
relationship between stress and strain in compression and to 
haviewiG ~eolstle Strength. she problem .is, analyzed.for .a 
prismatic member with a solid rectangular cross section with 


dimensions as shown in Figure l. 


i es 


Meiners in brackets refer to literature references listed 
in Appendix Tit. 
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Ste Loading Conditions 

Loading conditions are shown in Fig. 2. The compres- 
sive force P is a lineload which is evenly distributed 
along the width b of the member. The force is applied at 
an eccentricity equal to or greater than the kern-eccentricity 
Ca, oe Le Line or action of force P 1s parallel to 
the axis of the undeflected m2mber and the force is applied 
through hinges at both ends of the member. Thus the ends 
of the member are not restrained from rotation. The member 
is aiso not restrained from shortening in the direction of 


che Compressive force. 


4. EBqurlibraum 
Po EQualvbrium at -auCross Section 
hquUilabrium atsdny Cress Section.requires that. the 
resultant reaction force be equal to and collinear with the 
app Lied “force -P. ° The ‘Corresponding stress distribution 
Bet eaecLoss Section is OlLustratved in Pigs (3 for -two Gases. 
Figure 3(a) illustrates the case where P acts at an 
Seecitt Chey Cqudal Sto kermfeccentricity, te = t/6) "Figure 
3(b) illustrates the case where e > t/6. 
in the case allustrated in Figure 3(a) the compressive 
stress (co) at one “face (the tension rece oe the cross sec- 


tion is zero. At the other (compression) face, the maximum 


ee i ca 


—The term "tension face" is used here to identify the 
face that is subjected to zero stress and develops 
tension cracks, even though no tensile stress can exist 
in the member. 
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compressive stress (o ) occurs: 
o 


In Figure 3(b) a tensile crack appears at the tension 
Side of the cross section, since the material has no tensile 
strength. The uncracked part of the cross section has a 
triangular stress distribution similar to that shown in 
Fapure S(a), where o = 0 at the origin of the crack. 

The maximum stress at the compression face of the cross 


section is: 


2P (eq. 1) 


where u is the distance eS ene action tine of force. P 
ane the compression face :-of the cross section. 

the uncracked depth of the section in figure 3({b) is 
fa. oie ene depth of the crack is t-3Su. The expression in 
Pas Pets Valid tor all cases where e > t/6. This includes 


the case in figure 3(a), where u = t/3, and: 
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4:2 “Equilibrium*of theyStructural ‘Member 

The stress distribution within the structural member, 
Wienuokee Picvssapp ited tat ral ie¢ecentricity e¢ »>.t/6,; 1s 
illustrated in Fig. 4. The rectangle outlined by broken 
lines shows the undeflected shape of the member. The 
deflected shape is shown by the heavy outline. The shaded 
area within the deflected member shows the uncracked zone 
which supports the load. The stress distribution at one 
particular cross section is shown by the heavily -shaded 
triangle. 

The distance between the compression face of the mem- 
ber and the line of action of force P, u, varies along the 
height of the member, because of member deflection. The 
maximum distance, U,», occurs at the two member ends. The 
minimum distance, ul, occurs at midheight. The line of 
action of force P, which is shown in the figure by a dashed 
line, passes through the edge of the middle third of the 
compression zone. Since the stress at the compression face 
of each section, o> 1s inversely proportional to u, the 
maximum compressive stress in the member, o Occurs at 


max 
midheight: 


2P (eq. 2) 
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on The Differential Equation for the Deflection Curve 

Figure 5 shows the deflection curve of the compression 
“face of the member, together with the coordinate system used. 
Since the member is symmetrically loaded the origin was 
assumed at midheight. The x axis is parallel to the action 
line of P, and is tangential to the deflection curve at the 
Origin. At each point, u. = nena ae and, at x = h/2Z, 

a. thd 5 U,-: 

Figure 6 shows a small element of the deflected member 
between the compression face and the boundary of the uncracked 
zone, where o = 0 at all times. The element itself is shown 
by the solid outline. The average depth of the element is 
30... the initial undeflected’ length is d2. 

The curvature of the compression face is caused by the 
shortening of the face under compressive stress relative to 
the boundary of the uncracked zone, which does not change in 
length since o = 0 on the boundary. 

‘As d’+0, the following expression can be written for 
@, the change in the slope of the deflection curve over the 


length of the element: 


where € is the average strain over length df. 
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but a Feo eae (see eq. 1) 
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where o is the average stress at the compression face over 
° 


leneth al and 9 15 Young's modulus of elasticity, 


thus d = Zu 7 Sep. ae 


The. length of the compression face of the deflected 
element, d% (1 - €), can be expressed in terms of the radius 


6f curvature p: 


poo d2(1 = €) 
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The following equation can therefore be written for 


d*y/dx*: 
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This derivation is envisioned to apply to materials 
which develop relatively small strain (e < 0.005). The 
quantity 1 - e€ will therefore approximately equal unity. 
The maximum deflection developed will also be small rela- 
tive to the height of the member, and the slope of the 
deflection curve dy/dx, will be much smaller than 1 even 
at the point of maximum slope. The following approximation 
will therefore cause errors which are of lower-order 


magnitude and do not appreciably affect the result: 


3 
[1 = @y/ax)?y/2 
l-e i 


A similar assumption is traditionally made in derivations 
Ge G€lastic deflections and critical loads for conventional 
structural members. 


thus the expression for d*y/dx? reduces to: 


Since the value =r is constant for any particular member, 
the constant: 


ieee. 3 , : . 
Ky = ofp is substituted into the equation. 
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The following differential equation can therefore be written: 


el 
. (eq. 3) 
dx (uw +y) 


’ ay 4 ae 
where: Ky OED constant 


The following three boundary conditions apply to this equa- 


tron (refer to,.Fi9 S'S )* 


at x = 0, y = 0 fe 
ate: = 0., x 0 C2) 
ap k= A/2y cy = U, - U (3) 
G', Solution to they Differential -—Equation for the Deflection 
Curve 


Al €xacrysolutazon for qa 5S hasbeen derived:» Detailed 
steps of that derivation are presented in Appendix II. The 


following equation has been derived for the deflection curve: 
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ie DOUWdary,calgdition (3) :5 at b/Z; y= u, - u 


° 


is substituted into this equation, the following relationship 
between P, u,, and u, can be derived (See Appendix II): 
— = 0.40528 a| Ji-a+ta ln = + [2)] ° 
he : a ot 
where: 
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The term Bain is an "equivalent" critical load, computed from an 


requivalent" 1, which is based on an "equivalent" 


thickness of 3u,. 
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a as becomes the Euler load = when the section is 


loaded at the edge of the kern (3u,= t). 


° 
a = — 1s a measure of the magnitude of the maximum 
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deflection. 
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hs Application of the Solution to Derive Critical Load, 
| Maximum Deflection and Maximum Stress 


Figure 7 shows a graphical presentation of equation 5S. 
The plot is non dimensional and can be used to solve any 
particular numerical problem. 

Values of P/ Pie are plotted along the ordinate. On the 
abscissa, values of u /uy are plotted. It can be seen that 
each value of Pg? 3 corresponds to two values of ui /u,, 
except for the highest value of ieee that can be maintained 
by the member. The plot of the solution is interpreted as 
follows: | 

The value u/u, = 1 corresponds to an undeflected 
section. Thus this corresponds to no axial load or OPP ee = 0. 
Obviously, no axial load could be maintained if u. vanishes 
or u/uy = 0, since the mid section would be cracked over the 
full width, t. Thus the value Bye es = 0 also corresponds to 
u /uy = 0. 

When P is applied to an undeflected member and increased, 
deflections would increase correspondingly and the value of 
u /uy would decrease. This corresponds to the portion of 
the curve between u /u, = 1 and u /u, = 0.625, which repre- 
sents a condition of stable equilibrium. At the value of 
u /u, = 0.625 the maximum load is applied that the member 
can support, provided compressive strength was not exceeded 


at a lower value of ag Any further increase of P beyond 


"aa 


this point would cause the member to collapse. 
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Thus the important conclusion is reached, that elastic 


instability will occur at the critical load of: 


R = 0.285 Pia 


or Pax = 0.64 aes (eq. 6) 


The values of P/E ae. corresponding to values of u /u, in 
the range: 0.625 > u /uy, > 0 have little physical signifi- 
cance. They correspond to a second deflection that can be 
maintained by loads smaller than the critical load in a state 
of unstable equilibrium, 

Maximum stresses in the member, corresponding to any 
load smaller than the critical load ane can also be 
computed using the plot in Fig. 7 to compute ui» and eq. 2 


to compute maximum stress: 


eee: 
max 3bu 
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8. Summary of aieieeetca Analysis 

For the prismatic member shown in Fig. 1, made of 
material with no tensile strength, and subjected to the 
loading condition shown in Fig. 2. The following values 
can be computed from an exact solution to the differential 


equation for the deflection curve: 


<bne Crit cal load: 


Z 3 
T Ebu, 


geen = 0.64 am te 2 (eq. 6) 
* The maximum compressive stress: 


fe} = 2P (e 2) 
max bu cs 
° 


where u- can be determined from Fig. 7. 
° 


* The maximum deflection can be derived from Fig. 7. 
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9, . Application of the Theory to Predict the Results of 


Full-Scale Tests on Masonry Walls 


9.1 Source of Test Results and Test Specimens 

Only few test data are available which fit the boundary 
conditions and material properties assumed in this paper. 
The test results chosen for comparison were published by the 
Structural Clay Products Institute (SCPI) [2], and represent 
39 tests on full-scale brick walls of various slenderness 
ratios (h/t), loaded as shown in Fig. 2, Loads were applied 
at two eccentricities: e = t/6, and e = t/3. The walls were 
solid walls of single- and double-wythe construction, built 
in running bond, Slenderness ratios varied from h/t = 6.6 
to h/t = 46.1. Mortar was ASTM type S [1] and the compres- 
sive strength of the masonry oF) varied from 4,160 psi to 


5,100 psi. 


9.2 Comparison of the Assumed Material Properties with the 
Properties of the Test Specimens. 
While brick masonry has material properties similar to 
those assumed in the solution presented in this paper, several 
differences between the properties of the test specimens and 


the assumed idealized properties should be noted: 
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(1) The masonry used in the specimens tested should be 
expected to develop tensile strength of an order of 
Mapnitude of. 2-3 percent of its compressive strength. 
This property would somewhat change the equilibrium 
conditions illustrated in Fig. 3, which were assumed in 
the derivation. The change should result in somewhat 
greater-than-predicted capacities. The discrepancy 
would be greatest where stability failure occurs at 
‘relatively low pees: levels (large slenderness ratios 
and a eee 

£2Z),\..The stress-strain curve for brick is not exactly 
linear. The tangent modulus of elasticity at failure 
tends to be about 70 percent of the initial, tangent 
modulus of elasticity [3]. If deflections are predicted 
on the basis of the modulus of elasticity at low 

stress levels, deflections at high stress levels would 
probably be greater than the predicted deflections. 
This change in stiffness also slightly modifies the 
stress distribution shown in Fig. 3. These effects 
will increase with increasing stress levels at failure. 
(Small slenderness ratios and eccentricity). 

(3) The brick units themselves have greater strength 
and stiffness than the mortar beds connecting the units. 
This discontinuity results in a stress distribution 


which is much more complex than the idealized stress 
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distribution assumed in the solution. The net effect 
of these discontinuities on strength and stiffness 


has not been evaluated, 


9.3 Equations for Strength Calculations 

The test results presented by SCPI were obtained from 
a variety of brick walls, and for the sake of comparison, 
individual tests were normalized by dividing the vertical 
load carried by each wall by ae the load capacity of a 
cross section subjected to axial compression: io = f Bt. 
The value of fo? the compressive strength under axial load, 
was in each case derived from the average strength of short 
brick prisms which were constructed from similar brick and 
mortar at the time the corresponding test specimens were 
constructed and which were cured under similar conditions. 
It has been observed for brick masonry [3] that compressive 
strength in flexure exceeds the compressive strength in 
axial compression and that the strength increases with 
increasing strain gradients. In this case it has been deter- 
mined that cross sectional capacity of walls loaded at t/6 
and t/3 eccentricities can be approximately computed by using 
the stress distribution in Fig. 3 and assuming a compressive 


strength in flexure equal to 1.6 times the strength in axial 


compression: af. = 1.6f,. 
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It was also determined, that for the walls tested, the 
Moaudius, Ofpelasticaty is predicted with sufficient accuracy 
by assuming the fixed ratio of: file = ke Zo oe Mees 
The following equations can therefore be derived: 


Cross Sectional Capacity: 


ate. =Ft/oe- Pe = 068-72 
u fe) 


© 
ct 
o 

1) 


t/se. <P 0.4 P 
| O 
Wall Capacity: 
Walls may either fail when the flexural compressive 
strength is exceeded or when elastic instability occurs. 


For the case of a compression failure, wall capacity may 


be computed’ from ‘eq. 2: 


ee aty o>. P= 0.8 P. 
Uu fe) u, 
P u 
and — = 0.8 Bes 
Pp u 
Oo 1 
wr 
ae -G8= 0/33 P.=- 0.4 Pp. —— 
u fe) Uy 
a sae 
fe) 1 


u 
where — can be determined from Fig. 7. 


For a stability failure the following equations can be 


derived from eq. 6: 
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E = (£1 /1.215) 
and ps = fs bt 
Z. 2 30 
3 a it 
hs 192 Po( F) ( 5 
cr _ 192 (<2) 
Bern eje ns 
for e = t/6; 3u, = 
and er = 22 
ae SG xe 
Fe (h/t) 
for € = t/3; 3u, = 
and — per 
= are 
Pe (h/t) 


The following general 


a praction of Ps therefore 


at e = t/6: 
a i 
‘eee Oe 
fe) if 
gee yh 192 
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a) 
whichever is smaller 


at e = t/3: 
"u = 0.4 wo 
Po a 
sd Pu 24 
Pp, (h/t)? 


whichever is smaller 
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equations for wall capacity as 


apply: 


7) 


(eq. 


(eq. 8) 
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9.4 Comparison of Theory with Test Results 

Curves based on equations 7 and 8 are plotted in Fig. 
-8. Note that for the t/6 eccentricity compression failure 
occurs at slenderness ratios smaller than 22 (h/t < 22). 
For the t/3 eccentricity compression failure occurs at 
slenderness ratios smaller than 11 (h/t < 11). 

injhze. /oethe (theoretical curves are compared with 
test results. Test results at t/6 eccentricity are shown 
by the triangles and should be compared with the correspond- 
ing curve. Test results at the t/3 eccentricity are shown 
py ene Circles. 

himectcorar tie verena Of the test results is correctly 
predicted by theory and the magnitude of the load capacity 
teapproximately predicted. For the t/6 eccentricity capa- 
city was overestimated for h/t = 14 by about 10%. At 
h/t = 46 the capacity was underestimated by a substantial 
mate. probably because at this low stress level, the 
Petree OL the tensile strength is significant. For all 
other slenderness ratio the average capacity was approxi- 
mately predicted. For the t/3 eccentricity the capacity 
also tended to be underestimated at very low failure-stress 
levels, however, in the range of stability failure, the 
theoretical curve tends to approach the lower limit of the 


measured capacities. 
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9.5 Conclusion From the Comparison of the Theory with 
fest Results 
The theory correctly predicted the trend of the test 
results and approximately predicted the magnitude of mea- 
sured capacities. Since equations 7 and 8 are simple and 
relatively easy to use, their practical application should 


be seriously considered. 
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Append2x,1,.- Notation 


af” flexural compressive strength of masonry 
b width of member 

e load eccentricity 

dg incremental length 

E Young's modulus of elasticity 


fi, compressive strength of masonry under axial load 


h height of member 
Ky constant (Sep) : 
I moment of inertia = 
3 
ie equivalent moment of inertia — ) 
compressive force applied to member 
2 
crrercal, (Euler) load one 7) 
c Re 
T EI, 
- equivalent. critical ‘lee Cae 
a Critical load of member 
“ cross sectional axial-load capacity 
- load capacity under eccentric vertical compressive load 
£ thickness of member 
u distance between line of action of compressive load 
and compression face of member 
u distance between line of action of compressive load and 
° 
compression face of member at midheight 
Uj distance between line of action of compressive load and 


compression face of member at member support 
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dimensionless parameter = 


iL 
strain 
change in slope of deflection curve 
radius of curvature 
compressive stress 


maximum compressive stress in cross section 


— maximum compressive stress in member 
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APPENDIX II 


DERIVATION OF SOLUTION 


ae Ky ee 3 


Boundary Conditions: 


BE Ux = «0, y = 0 


at x = 0, X= 0 
x 
h 
at XK = 95> S darages: Damae 
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1 i 
Q = K fo” + u. luwas) -+ C,, ae ae inaju, + C., 
Co = -+° In Ju 
Ko o 
£ Vy + Vu, ty 
x -+E 5 ise) +. a, inf ae )| 


Substituting Ko 
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ie E Oo a/Uo 


Solving for u,: 
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Figure 1. Geometry and Dimensions of Member 
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Figure 2. 


Loading Conditions 
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Figure 3. Equilibrium at a Cross Section 
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Figure 5. Deflection Curve 
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Figure 6. Curvature 
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Figure 7. Graphical Presentation of Solution 
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Figure 8, Theoretical Strength of the Tested 
Walls as a Function of the Slenderness Ratio 
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Figure 9, Comparison of the Theory with 
Test Results 
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